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We investigate the quasi two-dimensional Taylor-Couette system in the regime where the radius 
ratio is close to unity - a transitional regime between three and two dimensions. By systematically 
increasing the Reynolds number we observe a number of standard transitions, such as one from the 
classical Taylor vortex flow (TVF) to wavy vortex flow (WVF), as well as the transition to fully 
developed turbulence. Prior to the onset of turbulence we observe intermittent burst patterns of 
localized turbulent patches, confirming the experimentally observed pattern of very short wavelength 
bursts (VSWBs). A striking finding is that, for Reynolds number larger than the onset of VSWBs, 
a new type of intermittently bursting behaviors emerge: burst patterns of azimuthally closed rings 
of various orders. We call them ring-burst patterns, which surround the cylinder completely but 
remain localized and separated by non-turbulent mostly wavy structures in the axial direction. We 
use a number of quantitative measures, including the cross-flow energy, to characterize the ring-burst 
patterns and to distinguish them from the background flow. The ring-burst patterns are interesting 
because it does not occur in either three- or two-dimensional Taylor-Couette flow: it occurs only 
in the transition, quasi two-dimensional regime of the system, a regime that is less studied but 
certainly deserves further attention so as to obtain deeper insights into turbulence. 

PACS numbers: 47.20.Ky, 47.32.cf, 47.54.-r 


I. INTRODUCTION 

Characteristics of turbulence in three- and two- 
dimensional flows are typically quite distinct. For ex¬ 
ample, in three-dimension flows, the energy spectrum 
of fully-developed turbulence obeys the well-known Kol¬ 
mogorov 1941 scaling law [I] of & -5 / 3 , while in two 
dimensions the scaling m m is k 3 . Although two- 
dimensional flow systems offer great advantages from the 
standpoint of computation and mathematical analysis as 
compared with three-dimensional flows, the former are 
not merely a kind of toy model of turbulence. In fact, 
two-dimensional turbulence is relevant to the dynamics 
of oceanic currents, origin of the ozone hole through mix¬ 
ing of chemical species in the polar stratosphere, the ex¬ 
istence of polar vortex, strong eddy motions such as trop¬ 
ical cyclones, and other large-scale motions of planetary 
atmospheres HE]. 

Turbulence is arguably one of the most difficult prob¬ 
lems in science and engineering, and the vast literature 
on turbulence was mostly focused on three or two di- 


*Electronic address: Sebastian.altmeyer@ist.ac.at 
^Electronic address: yhdo@knu.ac.kr 


mensions [6]. To gain new insights into turbulence, it 
is of interest to study the transitional regime between 
three and two dimensions. In such a regime, properties 
of both three- and two-dimensional flows are relevant, 
and one naturally wonders whether there are any unex¬ 
pected features associated with, for instance, transition 
to turbulence 0 0. The purpose of this paper is to 
report a new phenomenon in a paradigmatic quasi-two- 
dimensional system: the Taylor-Couette flow [9j with gap 
between the inner and outer cylinders so narrow that the 
system is neither completely three-dimensional nor ex¬ 
actly two-dimensional. In fact, this regime has not been 
investigated systematically previously. The new type of 
intermittent dynamics occurs en route to turbulence as 
the Reynolds number is increased. 

The Taylor-Couette system, a flow between two con¬ 
centric rotating cylinders, has been a paradigm in the 
study of complex dynamical behaviors of fluid flows, es¬ 
pecially turbulence fTOUTB] . The flow system can exhibit 
a large variety number of ordered and disordered behav¬ 
iors in different parameter regimes. For cylinders of rea¬ 
sonable length, the effective dimensionality of the system 
is determined by a single parameter - the ratio between 
the radii of the two cylinders. If the ratio is markedly 
less than unity, the flow is three dimensional. As the ra¬ 
tio approaches unity, the flow becomes two dimensional. 
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Most previous studies focused on the setting where the 
radius ratio is below, say about 0.9, the so-called wide- 
gap regime El, or when the ratio approaches unity so 
that the geometry is locally planar, resulting in an effec¬ 
tively two-dimensional Couette flow m- Here, we are 
interested in the narrow-gap case, where the radius ra¬ 
tio is close to unity but still deviates from it so that the 
flow is quasi two-dimensional. To be concrete, we fix the 
radius ratio to be 0.99 and, without loss of generality, 
restrict our study to the case where the outer cylinder 
is stationary. In fact, regardless of whether the outer 
cylinder is rotational or stationary, transition to turbu¬ 
lence can occur with increasing Reynolds number. In 
particular, for systems of counter-rotating cylinders, an 
early work m showed that transition to turbulence can 
be sudden as the Reynolds number is increased through 
a critical point. For fixed outer cylinder, the transition 
from laminar flow to turbulence can occur through a se¬ 
quence of instabilities of distinct nature nn [ 20 ]. 

For Taylor-Couette system of counter-rotating cylin¬ 
ders, spatially isolated flow patterns, the so-called lo¬ 
calized patches, can emerge through the whole fluid do¬ 
main and decay mm- Depending on the parameters, 
the localized patches can be laminar or more complex 
patterns such as inter-penetrating spirals. In the wide- 
gap regime (three-dimensional flow), numerical simula¬ 
tions [22J revealed the existence of so-called Gortler vor¬ 
tices [23] . small scale azimuthal vortices that can cause 
streaky structures and form herringbone-like patterns 
near the wall. Localized turbulent behaviors can arise 
when the Gortler vortices concentrate and grow at the 
outflow boundaries of the Taylor vortex cell [22^. 

For narrow gap (quasi two-dimensional) flows, there 
was experimental evidence of the phenomenon of very 
short wavelength bursts (VSWBs) [14]. One contribu¬ 
tion of our work is clear computational demonstration 
of VSWBs. Remarkably, we uncover a class of solutions 
in quasi two-dimensional Taylor-Couette flow en route 
to turbulence. These are localized, irregular, intermit¬ 
tently bursting, azimuthally closed patterns that man¬ 
ifest themselves as various rings located along the ax¬ 
ial direction. For convenience, we refer to the states as 
“ ring-bursts The number of distinct rings can vary de¬ 
pending on the parameter setting but their extents in 
the axial direction are similar. The ring bursts can occur 
on some background flow that is not necessarily regu¬ 
lar. For example, in typical settings the background can 
be wavy vortex flows (WVFs) with relatively high az¬ 
imuthal wave numbers. Because of the coexistence of 
complex flow patterns, to single out ring bursts is chal¬ 
lenging, a task that we accomplish by developing an ef¬ 
fective mode separation method based on the cross-flow 
energy. We also find that ring bursts are precursors to 


turbulence, signifying a new route to turbulence uniquely 
for quasi two-dimensional Taylor-Couette flows. To our 
knowledge, there was no prior report of ring bursts pat¬ 
terns or likes. This is mainly due to the fact that the 
quasi two-dimensional regime is a less explored territory 
in the giant landscape of turbulence research. It would be 
interesting to identify precursors to turbulence in quasi 
two-dimensional flow systems in general. 

In Sec. [n] we outline our numerical method and 
describe a number of regular states in quasi two- 
dimensional Taylor-Couette flows. In Sec. |III| we present 
our main results: numerical confirmation of experimen¬ 
tally observed VSWBs and more importantly, identifi¬ 
cation and quantitative confirmation of intermittent ring 
bursts as precursors to turbulence. In Sec. E| we present 
conclusions and discussions. 


II. NUMERICAL METHOD AND BASIC 
DYNAMICAL STATES OF QUASI 
TWO-DIMENSIONAL TAYLOR-COUETTE 
FLOW 

A. Numerical method 

The Taylor-Couette system consists of two indepen¬ 
dently rotating cylinders of finite length L and a fluid 
confined in the annular gap between the two cylinders. 
We consider the setting in which the inner cylinder of ra¬ 
dius Ri rotates at angular speed Q and the outer cylinder 
of radius R 0 is stationary. The end-walls enclosing the 
annulus in the axial direction are stationary and the fluid 
in the annulus is assumed to be Newtonian, isothermal 
and incompressible with kinematic viscosity v. Using the 
gap width d = R 0 — Ri as the length scale and the ra¬ 
dial diffusion time (DT) d 2 jv as the time scale, the non- 
dimensionalized Navier-Stokes and continuity equations 
are 


d t u + (u • V)u = —Vp + V 2 u, V • u a» 0, (1) 

where u = (u ri uq^ u z ) is the flow velocity field in cylin¬ 
drical coordinates (r, 0, z), the corresponding vorticity is 
given by V x u = (£, 77, £), and r is the normalized radius 
for fluid in the gap (0 < r < 1). The three relevant pa¬ 
rameters are the Reynolds number Re = the ra¬ 

dius ratio Ri/R 0 = 0.99, and the aspect ratio T = L/d = 
44. The boundary conditions on the cylindrical surfaces 
are of the non-slip type, with u(ri,0,z,t) = (0,i?e, 0), 
u(r o ,0, z,t) = (0,0,0), where the non-dimensionalized 
inner and outer radii are r* = Ri/d and r Q = R 0 /d , 
respectively. The boundary conditions in the axial direc¬ 
tion are u(r, 0 , T0.5T, t) = (0, 0, 0). 

We solve Eq. 0 by using the standard second-order 
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time-splitting method with consistent boundary condi¬ 
tions for the pressure [24]. Spatial discretization is done 
via a Galerkin-Fourier expansion in 6 and Chebyshev col¬ 
location in r and z. The idealized boundary conditions 
are discontinuous at the junctions where the stationary 
end-walls meet the rotating inner cylinder. In experi¬ 
ments there are small but finite gaps at these junctions 
where the azimuthal velocity is adjusted to zero. To 
achieve accuracy associated with the spectral method, a 
regularization of the discontinuous boundary conditions 
is implemented, which is of the form 

u e {r,6,±Q.5T,t) = i?e{exp([rj - r]/e) + exp([r - r 0 ]/e)}, 

(2) 

where e is a small parameter characterizing the physical 
gaps. We use e = 6 x 10 -3 . Our numerical method 
was previously developed to study the end-wall effects in 
the Taylor-Couette system with co- and counter-rotating 
cylinders [25j 26 . In the present work we use up to n r = 
50 and n z = 500 Chebyshev modes in the radial and 
axial directions, respectively, and up to ng = 100 Fourier 
modes in the azimuthal direction. The time step is chosen 
to be St = 10 -6 . 

B. Qualitative description of basic dynamical states 

1. Low order instabilities 

In the wide-gap, three-dimensional Taylor-Couette sys¬ 
tem, various flow patterns and their bifurcation behav¬ 
iors are relatively well understood mm In our quasi 
two-dimensional setting, a number of known low-order, 
non-turbulent instabilities still occur, which constitute 
the background flow upon which a new type of ring burst 
structures emerge. Here we describe these low-order in¬ 
stabilities. To visualize and distinguish qualitatively dif¬ 
ferent flow patterns, we use the contour plots of the az¬ 
imuthal vorticity component g. 

Primary instability - Taylor vortex flow. As the 
Reynolds number Re is increased, the basic state, circu¬ 
lar Couette flow (CCF) becomes unstable and is replaced 
by the classical Taylor vortex flow (TVF) that consists 
of toroidally closed vortices. TVFs with some increasing 
numbers of vortices appear gradually over a large range 
of Re, starting from a single Taylor vortex cell generated 
by the mechanism of Ekman pumping near the end walls. 
This initial cell can appear either near the top or the bot¬ 
tom lid. In a perfectly circular geometry, initial Taylor 
vortex cells can occur simultaneously at both lids, but 
this is less likely in realistic systems due to the inevitable 
imperfections in the system. The onset of initial TVF 
cells was experimentally found for Re about [14] 358. We 
find, numerically, that the onset value is about 356. As 


Re is increased from this value, additional vortices ap¬ 
pear, which enter the bulk from near the lids until they 
finally fill the whole annulus for Re value about 435 (ex¬ 
perimentally the value is about [141 437). For example, 
we observe a TVF with 22 pairs of vortices within the an¬ 
nulus, with the characterizing wavenumber of k = 3.427, 
and an additional pair of Ekman boundary layer vortices 
near the top and bottom lids. 

Secondary instabilities - wavy vortex flow. Further in¬ 
creasing Re, the TVF becomes unstable and is replaced 
by wavy vortex flow (WVF), a kind of centrifugal insta¬ 
bility that appears via a supercritical Hopf bifurcation, 
which exhibits a wavy-like modulation but the vortices 
remain azimuthally closed. We find that the onset of 
WVF occurs for Re ze 473 (experimentally the value is 
about |14] 475). The WVF, as shown in Fig. [lja), is as¬ 
sociated with relatively high values of the azimuthal wave 
number m. In general, WVF states consisting of many 
azimuthal modes can occur. For example, the WVF pat¬ 
tern in Fig[lja) contains m = 39 modes, which we name 
as WVF 39 . We find that this particular WVF state is in 
fact a global background flow state in a wide range of Re 
values. 

The TVF and WVF patterns occur regardless of 
whether the flow is three dimensional or quasi two- 
dimensional, although the critical values of Re for the 
onsets of these flow patterns depend on the system de¬ 
tails. 


2. High-order instabilities 

Very short wavelength bursts and localized patches. 
For the commonly studied [13], wide-gap, three- 
dimensional Taylor-Couette system, e.g., of radius ra¬ 
tio 0.883, the typical sequence of solutions with increas¬ 
ing Re values is as follows. As the stationary TVF be¬ 
comes unstable, time dependent WVF arises, followed 
by the so-called modulated WVF m, eventually lead¬ 
ing to turbulent behaviors. However, in the narrow-gap, 
quasi two-dimensional system, we have not observed the 
global transition from WVF to modulated WVF. Instead, 
we find that very short wavelength bursts (VSWBs) ap¬ 
pear directly after the onset of WVF without any other 
types of intermediate solutions. Onset of VSWBs is 
Re ~ 483, which agrees with the experimentally observed 
onset value m- 

As Re is increased, we observe localized patches (LPs) 
superposed on WVF. The pattern within LPs can be ei¬ 
ther wavy-like or turbulent, depending on the Re value 
(lower and high values for the former and the latter, re¬ 
spectively). In particular, LPs with wavy-like interior 
behaviors emerge, and for Re > 580 LPs with turbu- 
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(b) Ring bursts (n=3) for Re=545 





(c) VSWB for Re=800 




FIG. 1: (Color online) Distinct flow patterns in quasi 
two-dimensional Taylor-Couette system. Contours of 
the azimuthal vorticity component 77 with (left column, 
7 / E [—400,400]) and without (right column, 77(771 ^ 0) E 
[— 200 , 200 ]) the axisymmetric component on an unrolled 
cylindrical surface at the midgap (r = d/2) for different flows 
and values of Re. Red (dark gray) and yellow (light gray) col¬ 
ors correspond to positive and negative values, respectively, 
and the black line indicates the zero contour. 


lent behaviors arise, which can either grow to VSBWs 
or decay slowly. The number of LPs depends on Re as 
well. In general, the higher the value of Re the larger 
the number of LPs appears. The patches are randomly 
distributed over the whole bulk length and are always 
visible except for the region near the lids in which strong 
Ekman vortices arise. 

Ring-burst patterns. When Re is increased above 
about 540, we discover a new type of localized, inter¬ 
mittent burst solution, the ring bursts that coexist with 
VSWBs. While both types of solutions are localized, 


there are characteristic differences. For example, VSWBs 
appear randomly over the whole bulk fluid region with 
seemingly expanding behaviors in all directions, but ring 
bursts always remain localized in the axial direction. In 
fact, the bursts are generated from the localized turbu¬ 
lent patches that grow in the azimuthal direction as Re is 
increased. For sufficiently high values of i?e, the patches 
extend over one circumference, generating distinct ring 
bursts that are separated from the flow patterns in the 
rest of the bulk, as shown in Fig. [ljc) . The character¬ 
istics of the background flows in the regions surround¬ 
ing the ring bursts depend strongly on Re. They can 
range from wavy-like patterns (i.e., WVF 39 ) to inter¬ 
penetrating spirals due to the interactions among various 
azimuthal modes, as shown in Fig. [ljc). We observe ring 
burst patterns of different order n, as shown in Figs, [l] 
and [2|. For example, we observe n E {1,2,3}. All these 
states coexist but the probability for ring bursts with 
larger values of n increases with Re. Along the axial 
direction, the ring burst patterns can appear at any po¬ 
sition, except for the regions of Ekman vortices near the 
lids, because strong boundary layer vortices prevent the 
development of ring burst patterns. 

While there is no apparent order associated with the 
flow patterns within the ring bursts, the flow in the re¬ 
gions in between exhibit a clear signature of WVF 39 pat¬ 
tern. The ring bursts can spontaneously break up and 
disappear. Depending on the number of ring bursts, the 
transient time that it takes for the burst to decay into lo¬ 
calized patches can be relatively long. In addition, there 
can be transitions between patterns with distinct num¬ 
bers of ring bursts. For example, suppose there is a ring- 
burst pattern of order n = 3. If one ring disappears, 
a new ring-burst region can appear and grow. Transi¬ 
tions between patterns with either increasing or decreas¬ 
ing numbers of ring bursts have been observed. 

For Re > 850, we find that VSWBs fill the whole bulk 
fluid region, as shown in Fig. [TJc). In fact, for Re > 900, 
the whole fluid region is saturated with VSWBs. These 
numerical observations are in agreement with experimen¬ 
tal findings El- Note, however, even when VSWBs fill 
the interior of the bulk WVF-like patterns always appear 
near the axial end-walls due to the fact that Ekman vor¬ 
tices can stabilize the boundary layers against turbulent 
bursts. The size of the WVF-like Ekman regions depends 
on Re - for higher value of Re the smaller the WVF re¬ 
gions near the lids appear, as can be seen from Fig. [TJ 

Underlying flow patterns. The right panels in Fig. [l] 
show contours for the same pattern as those for the left 
panels but without the underlying axisymmetric contri¬ 
bution. The resulting “reduced” flow pattern of WVF 39 
in Fig. [ija) appears quite regular, indicating the exis¬ 
tence of the dominant wavy mode with m = 39 in the 
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FIG. 2: (Color online) Angular momentum and vortic- 
ity. Isosurfaces of rv (top row) and 77 (bottom row) for flows 
at different values of Re (isolevels shown for the top and bot¬ 
tom rows are rv = 80 and 77 = ±30, respectively). For clear 
visualization, here and in all following three-dimensional plots 
the radius ratio is scaled by the factor 100 . 


azimuthal direction. The black curves in the vertical 
direction represent the zero contour lines. At several 
azimuthal positions, these lines narrow, indicating the 
emergence of a second but weaker mode, e.g., m = 8 . 
For order 3 ring bursts, the flow patterns near the lids 
are somewhat modified due to the presence of higher m 


modes, as shown in Fig. [ljb). The flow patterns in the 
central region (including that containing the three ring 
bursts) exhibit a completely different behavior. In par¬ 
ticular, with respect to contours 77(771 7^ 0) the turbulent 
ring bursts and the separating WVF 39 pattern appear 
indistinguishable. This finding suggest a quite strong ax- 
isymmetric dominance of the flow even when ring bursts 
arise, confirming that the bursting structure is ordered 
in azimuthally closed rings around the cylinder with the 
same m = 0 axisymmetry. For turbulent flows or VSWBs 
in Fig. 0c), none of the patterns has such an ordered 
structure. 


III. CHARACTERIZATION OF RING BURSTS 

A. Angular momentum, azimuthal vorticity, and 
modal kinetic energy 


Figure |2] shows the isosurfaces of the angular momen¬ 
tum rv (top row) and azimuthal vorticity 77 (bottom row) 
for representative flow patterns at different values of Re. 
The dominant axisymmetric contribution and high az¬ 
imuthal wavenumber associated with WVF 39 , as shown 
in Fig. [2p), serve as the background pattern for ring 
bursts and VSWBs. The three turbulent, azimuthally 
closed burst regions associated with the order-3 ring 
burst pattern is distinctly visible, as shown in Fig. |2|b). 
Within each burst region, both rv and 77 appear random 
but the (background) flow patterns in between the burst¬ 
ing regions are remnant of the WVF 39 pattern. With 
increasing Re higher m modes emerge, but the separa¬ 
tion between burst and non-burst regions persists. For 
high values of Re, VSWBs arise, as shown in Fig. |2|d). 
In this case, isosurface plots of rv and 77 exhibit random 
flow patterns except for the regions near the lids due to 
the Ekman boundary layers. 

To better visualize and illustrate the similarity and 
differences between ring burst and VSWB patterns, we 
present in Fig. [3] the behaviors of the angular momentum 
(top panels) and azimuthal vorticity (bottom panels) for 
a segment of the bulk flow. The length of the segment 
is T/IO. We see that, for the ring burst pattern (left 
column), there is a dominant axisymmetric component, 
but this is lost for VSWB. While the entire structure of 
the ring-burst pattern possesses an axisymmetry due to 
its azimuthal dominance (m = 0 ), no apparent symmetry 
exists in the interior of the burst regions. There are thus 
two distinct spatial scales associated with ring bursts: a 
large scale determined by the axisymmetry and a small 
scale present in the interior of the bursting regions. 












6 


(a) Ring bursts of order 3 (b) VSWB 

Re=545 Re=800 






FIG. 3: (Color online) Three-dimensional views of angu¬ 
lar momentum and vorticity. Magnitude of the angular 
momentum (top panels) and azimuthal vorticity (bottom pan¬ 
els) for a flow segment of vertical length T/IO for ring bursts 
(left column) and VSWB (right column). The segment is 
taken from the region containing the middle ring burst in 
Fig- Ip). 
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FIG. 4: (Color online) Behaviors of the kinetic energy. 

(a) Time-averaged (« 50 diffusion time) axisymmetric en¬ 
ergy E o versus Re and (b) time series of E 0 for a number 
of representative values of Re. The time intervals in which 
distinct ring-burst patterns last are indicated. For exam¬ 
ple, for Re = 560, ring burst pattern of order 3 is present 
for 2.5 < t < 6.9, and that of order 2 can be found for 
9.2 < t < 15.8 and 42.6 < t < 47.1. For Re m 545, order-2 
ring burst pattern appears for 0.5 < t < 4.5 and the order-3 
pattern is present for 40.3 < t < 44. However, the ring burst 
patterns are not reflected in Eo(t). 


We next examine the modal kinetic energy defined as 


_ r 27t nT nr 0 

E := T] E m = / / / u m u* m rdrdzd 6 >, 

™ Jo Jo Jn 


(3) 


where u m (u* m ) is the ra-th (complex conjugated) 
Fourier mode of the velocity field. Due to time depen¬ 
dence of the solutions it is necessary to calculate the time 
averaged kinetic energy E. In the energy characterization 
local quantities such as the radial velocity at mid-height 
and mid-gap, u r (d/2, 0 , T/ 2 , £), are relevant. 

Figures [4^a) andQb) show the time-averaged axisym¬ 
metric energy E o versus Re and its time evolutions for a 
number of representative Re values, respectively. We see 
that E o increases monotonously with i7e, regardless of 
the nature of the bulk flow patterns. Likewise the time 
evolution of the kinetic energy presents no clear indica¬ 
tion of ring burst patterns. For example, for Re = 545, 
in the time interval [40,44] the flow pattern should be 
order-3 ring burst, but the energy time evolution does 
not show any signature of this burst pattern. 

Figures [5ja) and §b) show the power spectral den¬ 
sity (PSD) of the radial velocity profile at the midgap 
for WVF and ring burst, respectively. We see that the 
PSD associated with the ring burst pattern (b) indicates 
the existence of significantly higher modes than the WVF 
pattern (a). In fact, the PSD of WVF 39 in (a) shows a 
strong peak at the frequency of about 19, which corre¬ 
sponds to the dominant azimuthal mode (m = 39). This 
peak is still present in (b) but it is much broadened, indi¬ 
cating WVF as the background flow pattern for the ring 
burst. Figure ^c) shows the scaled PSD curves for sev¬ 
eral flow patterns for different values of Re. We see that 
the PSD curves for the ring-burst patterns essentially col¬ 
lapse into one frequency band, but the PSD curve asso¬ 
ciated with WVF 39 lies slightly below those of ring burst 
patterns. There is relatively large difference for small 
frequencies but it becomes insignificant for higher fre¬ 
quencies. This is further support for the role of WVF 39 
in providing the skeleton structure for all ring burst pat¬ 
terns. 


B. Cross-flow energy 

From Fig. |4](b) , we see that the modal kinetic energy is 
not effective at distinguishing ring bursts from other basic 
flow patterns. A suitable and commonly used quantity 
in the study of fluid turbulence is the so-called cross-flow 
energy, where typically small values indicate a flow being 
laminar, large values corresponding to turbulence. 
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FIG. 6 : (Color online) Radial component of cross-flow 
energy. Spacetime plots of the radial component of the 
cross-flow energy, E c ^ r (rfl) = (u% + u 2 z )A( r )i averaged over 
the surfaces A of a concentric cylinder at radius r for (a) 
Re = 500, (b) Re = 545, and (c) Re = 560. Red (yellow) 
color indicates high (low) energy with contours defined as 
AE cf ' r = 5 x 10 5 . The maximum energy values in (a-c) are 
approximately 8.1xl0 5 , 8.5xl0 6 , and 1.2 x 10 6 , respectively. 
For better visualization the radial gap width is magnified by 
a factor of 500 (the same for Fig. [7] below). 


FIG. 5: (Color online) Power spectral density (PSD). 
PSD curves calculated from the radial flow component Ud /2 = 
u(d/ 2,0,r/2,t) for (a) WVF 39 (Re — 500) and (b) ring burst 
of order 3 for Re = 545. (c) PSD curves scaled by the re¬ 
spective Reynolds number for different flow patterns. As the 
number of ring bursts is increased, the local peaks in the PSD 
curves become more pronounced. 


1. Radial component of cross-flow energy 

The radial component of the cross-flow energy is given 
by [27] 

E cf,r ( r,t) = (u 2 r + u 2 z ) A(r) , (4) 

where (-) A(r) denotes averaging over the surface of a con¬ 
centric cylinder at radius r. The cross-flow energy com¬ 
ponent E c f ,r (r,t) measures the instantaneous energy as¬ 
sociated with the radial and axial velocity components at 
radial distance r. Figure [ 6 ] shows the spacetime plots of 
E c f’ r (r, t ) over the time period of 50 DT for three values 
of Re. We observe the temporal emergence and disap¬ 
pearance of various ring-burst patterns in the bulk. For 
example, for Re = 545, the order-3 ring burst pattern 
exists for 40 < t < 44, and the order-2 pattern appears 
for 0.5 < t < 4.5, as shown in Fig. [ 6 ^b). For Re = 560, 
an order-3 pattern appears for 2.5 < t < 7.0, and a sin¬ 


gle ring burst pattern (order 1 ) exists for 43 < t < 47, as 
shown in Fig. [ 6 jc). 

Plots of the cross-flow energy exhibit two features. 
Firstly, presence of the order -1 ring burst pattern is ac¬ 
companied by a significant increase in the radial cross- 
flow energy as compared with that associated with the 
background flow, e.g., WVF 39 , as indicated by the uni¬ 
form red (dark gray) regions. Secondly, the profile of 
E c f’ r (r,t) for any ring burst pattern is approximately 
symmetric with respect to the middle of the gap. It is 
thus plausible that the ring-burst patterns are similar to 
the so-called turbulent WVF state [ 22 , [281430] (see Sec. 


Ill C). The symmetry may be a consequence of the quasi 


two-dimensional nature of the flow. In all cases, the oc¬ 
currence of ring burst patterns enhances the magnitude 
of the cross-flow energy. 

Figure [7][a) shows a magnification of a segment of 
Fig. ^b) for t G [0,1.8], which contains the emergence 
and disappearance of order-3 ring burst for Re = 545, 
occurring at t & 0.45 and t ~ 1.35, respectively. Outside 
this time interval the flow is WVF 39 , which constitutes a 
nearly uniform background without any significant vari¬ 
ation in E c b r (r, t). The spatial distribution of E c f' r (r, t) 
exhibits a kind of symmetry about r = 0.5. Figure [7][b) 
shows the average cross-flow energy, (E c f ,r ) r , as a func¬ 
tion of time £, corresponding to the emergence and disap- 


























FIG. 7: (Color online) Magnified view of the radial com¬ 
ponent of the cross-flow energy, (a) Magnification of a 
segment of Fig. |6jb) for t £ [0,1.8] illustrating the emergence 
and disappearance of order-3 ring burst pattern for Re = 545. 
The contours are defined through A E cf,r = 5 * 10 5 , and the 
maximum energy value is about 8.288 x 10 6 . (b) Averaged ra¬ 
dial cross-flow energy ( E cf,r ) r versus t for order-3 ring burst 
for Re = 545. The nearly constant background flow is WVF39 
with (£ c/ ’ r (WVF 39 ))r « 2 x 1 (T 6 . 


pearance of the ring burst pattern. We observe a signifi¬ 
cant enhancement of the average radial energy over that 
of the background flow {{E c ^ r ) r ^ 2 x 10 -6 . In fact, in 
the time interval where the ring burst exists, the maxi¬ 
mum value of the average radial energy is about one order 
of magnitude larger than that of the background flow. In 
general, the maximum energy depends on the number n 
(order) of ring bursts in the annulus, where a larger value 
of n corresponds to larger value of the maximum energy. 

The emergence and development of any type of burst 
patterns with increasing Re can be conveniently charac¬ 
terized using the quantity /#, the percentage of the an¬ 
nulus showing bursts in the spacetime plot. Figure |8ja) 
shows /b versus Re, where we observe approximately 
a linear behavior for Re < 900, and the transition to 
bursts occurs for Re « 480, in agreement with the on¬ 
set of VSWB |T4]. Onset of ring bursts can be revealed 
by examining the maximum value of the average radial 
energy, (E c f ,r ) r , versus Re, as shown in Fig. ^b). Re¬ 
gardless of the order n of ring bursts, there is a linear in¬ 
crease in (E c f ,r ) r with Re, suggesting a type of forward 
bifurcation. Calculations of the flow amplitudes show 
a square-root type of scaling behavior with increasing 
parameter difference from the “critical” point, providing 
further support for the forward nature of the bifurcation. 
Due to strong localization of the ring bursts, the slope of 
the linear scaling behavior depends on the ring-burst or¬ 
der n. Nonetheless, the onset value Re c ~ 537 of ring 



FIG. 8: (Color online) Onset of burst patterns, (a) Burst 
fraction R as a function of Re. A linear fit indicates that 
onset of VSWBs occurs for Re « 480 (the experimental on¬ 
set m is Re ~ 483). (b) The maximum average radial energy 
(E cf ’ r )r [Eq. @] versus Re for ring bursts of different order 
(n £ {1,2,3}). The linear behavior suggests that ring bursts 
are result of a forward bifurcation, the onset of which occurs 
for Re c ~ 537. In this figure, the average energy of the un¬ 
derlaying background flow, ( E cf,r ) r , has been subtracted off. 

bursts does not depend on n. This is evidence that ring 
bursts emerge after VSWBs. 


2. Axial component of cross-flow energy 
The axial component of the cross-flow energy is 

E^' = (( U r)m V { u z)m)A(z)- (5) 

where A(z) stands for averaging over the radial and az¬ 
imuthal variables on the surface of a disc at a fixed axial 
position z. Figured shows the time-averaged value E m 
for the order-3 ring burst pattern for different azimuthal 
modes. The axisymmetric component E 0 ’ is dominant. 
Near the center of the ring burst region in the axial direc- 
tion, the values of E 0 ’ are smaller than those around the 
edges. However, the contributions from all higher modes, 
i.e., E m for m > 0, have similarly small magnitudes as 
compared with E°J ,Z , as shown in Fig. [^b). In fact, for 
higher m modes, the axial components of their cross-flow 
energies are randomly distributed over z and they are 
not indicators of any appreciable difference between the 
background flow and the ring burst pattern. 
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FIG. 9: (Color online) Axial component of cross-flow 
energy. Associated with the order-3 ring burst pattern for 
Re — 545, (a) time averaged axisymmetric energy component 
E o ’ and (b) energy contributions E m ’ from selected higher 
azimuthal m modes. In (b), additional axial average energy 
values are indicated. 

3. Mode separation of axial cross-flow energy 

To better characterize the ring-burst patterns in rela¬ 
tion to the background wavy-like and general burst pat¬ 
terns, we devise a method based on the idea of mode sep¬ 
aration. Since a burst pattern includes modes of higher 
azimuthal wave numbers, we can decompose the axial 
component of the cross-flow energy E^ z {zfl) into two 
distinct subcomponents: 

m c — 1 M 

E c J’*(z)+E c b f ’ z (z) = E%>*{z)+ E *#■*(*)> (6) 

m=1 m=m c 

where Eff' z (z) and E^' z (z) denote the axial components 
of the cross-flow energy associated with the background 
wavy-like and burst patterns, respectively, and m c is 
some cutoff mode number. Note that the axisymmet¬ 
ric component of the cross-flow energy is excluded be¬ 
cause it is significantly larger than all other components 
[Figs. [9^a,b)]. We choose the normalization factor to be 
the total cross-flow energy for all modes except m = 0: 

M 

e m z = E E °d ,z - (?) 

m= 1 

Figure [lQ|[a) shows the basic E w and the burst contri¬ 
bution Ei) at different axial positions versus the cut¬ 




FIG. 10: (Color online) Characterization of axial cross- 
flow energy based on mode separation, (a) Varia¬ 
tion of the axial cross-flow energy component E^ ,z (z) = 
E c b f ’ z (z) + E c J’*(z) versus the cutoff wavenumber m c for two 
axial positions in the burst region (z = 0.3r and 0.45r, cir¬ 
cles) and in the background region ( z = 0.1 5r and 0.35r, 
squares), (b) Spatial variations of E^ f,z {z) = E m ,z ( z ) 

and Eyf ,z (z) J2m=m c E ™ ,Z ( z ) with ^ for m c = 20 (circles) 
and 30 (squares). Regions above (below) each curve indi¬ 
cate Eb (E w ). The two horizontal dashed lines give the 80% 
threshold of the maximum values of [E w (m c == 20)] (lower) 
and [E w (m c — 30)] (upper), respectively. In the calculations 
the zero mode contribution is excluded. Due to normalization 
with Em the sum must be unity. 

off wavenumber m c . We observe that the curves for 
z positions within the burst region (circles, e.g, at 
z = 0.3T and0.45T, Figs, [l] and [9| are higher than 
those in the wavy-like background (squares, e.g., at 
z = 0.15r and 0.35r), where the former exhibit a rapid 
decrease in the energy to collapse with the latter for m c 
about 39. The variations of E 5 and E w along the annulus 
length for two different cutoff wavenumbers [m c = 20, 30, 
as indicated by the vertical lines in Fig.[Io|[a)] are shown 
in Fig. mb). Neglecting the differences in their mag¬ 
nitudes, the variations show qualitatively similar behav¬ 
iors. Using a cutoff level with 80% of the maximum of 
the background contribution [. E w (m c )\ , we find a good 
agreement with the visible energy thresholds between the 
background and burst regions (cf. Fig. [l]). We thus see 
that, through some proper choice of the cutoff mode num¬ 
ber, the burst and non-burst regions can be distinguished 
by examining the axial cross-flow energy variations asso¬ 
ciated with the cutoff mode. 
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FIG. 11: (Color online) Axial spacing of ring bursts, (a) 

SPSD of the radial velocity u r at the midgap for Re = 545, 
where k is the axial wave number, (b) Axial spacing of the 
burst regions for order-n ring bursts (n G {1,2, 3}), indicated 
by the horizontal lines. The average spacing is about 0.079, 
the axial spacing of two pairs of Taylor vortices constituting 
four vortex cells. 

C. Axial spacing and localization 

The axial ranges of distinct ring burst regions are ap¬ 
proximately identical. Figure [llja) shows the spatial 
power spectral density (SPSD) of the radial velocity u r 
along a line at the midgap in the center of the bulk for 
the background flow and three types of ring-burst pat¬ 
terns. We observe that the four curves coincide at the 
first sharp peak determined by the wave number associ¬ 
ated with the background flow WVF 39 that consists of 
24 vortex pairs in the axial direction. The correspond¬ 
ing axial wavelength and wave number are A ~ 1.667 
and k « 3.770, respectively. In addition, several broad¬ 
band peaks at higher wave numbers exist in the SPSD 
of the ring burst patterns, corresponding to a number of 
short wavelength bursts within the respective patterns. 
The broadband nature at higher wave numbers are in¬ 
dicative of dominance of small-scale burst patterns. Fig¬ 
ure [Tljb) shows the axial spacing associated with various 
ring burst patterns as a function of Re, where each hor¬ 
izontal dashed line represents the averaged axial spacing 
for a particular ring burst pattern. The axial spacing 
is apparently independent of the value of Re and of the 
order n of the ring-burst pattern. The typical value of 
axial expansion agrees well with the axial dimension of 


two-pair Taylor vortices that constitute four single vor¬ 
tex cells, which holds for all ring-burst flows that we have 
succeeded in uncovering. Analogous to the behavior of 
the cross-flow energy, this behavior is indicative of that 
of turbulent WVFs j22 1 i28l I3U] . We note that the size of 
only one pair of Taylor vortices (two cells) is too small to 
account for the observed range of axial expansion. This 
is consistent with the formation process of ring bursts. In 
particular, any localized turbulent patch, after its gener¬ 
ation, first expands in the axial direction (to four cells) 
before growing in the azimuthal direction. Whenever the 
burst region has expanded to a larger size in the axial 
direction, the closed ring structure is destroyed, leading 
to VSWBs. 


IV. CONCLUSIONS AND DISCUSSIONS 

This paper provides a comprehensive numerical study 
of the quasi two-dimensional Taylor-Couette systems of 
radius ratio close but not equal to unity, a regime that 
has not been studied previously. The relevant control or 
bifurcation parameter is the Reynolds number Re. For 
small Re values, TVF initially arises near one of the lids 
in a single cell, extends and finally fills the bulk inte¬ 
rior completely. As Re is increased the TVF loses its 
stability and WVF emerges through a supercritical Hopf 
bifurcation. WVFs with high azimuthal wavenumbers, 
e.g., m — 39, constitute a persistent background flow, on 
top of which more complex flow structures develop, such 
as VSWBs that has been experimentally observed M- 

The main result of this paper is the discovery of a new 
type of transient, intermittent state en route to turbu¬ 
lence with increasing Re: ring bursts. They emerge when 
localized turbulent patches grow and close azimuthally, 
signifying a higher order instability. The ring bursts oc¬ 
cur in various orders and the azimuthally closed burst 
regions possess axial expansion surrounded by the WVF 
background flow or more complex flows at high Re values. 
The ring burst patterns differ characteristically from the 
localized VSWB turbulent patches. The axial expansion 
of ring burst patterns of different orders correlates well 
with the size of the double-pair Taylor vortex structure, 
providing plausible reason that the patterns are similar 
to the turbulent WVF structure [22] 1281130] that usu¬ 
ally occurs through the whole annulus at higher Reynolds 
numbers. We develop a mode separation method based 
on decomposing the cross-flow energy to distinguish be¬ 
tween burst and non-burst patterns, where the former 
and latter are associated with higher and lower azimuthal 
modes, respectively. We also find that the radial cross- 
flow energy changes significantly in the presence of ring 
burst patterns. By examining the maximum value of 
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the cross-flow energy, we determine the onset of the 
ring-burst patterns at the critical Reynolds number of 
Re c ~ 537, which is larger than that for the onset of 
VSWBs [14] (about 482). For ring burst patterns of dif¬ 
ferent orders, their expansions in the axial direction are 
nearly identical, and they tend to shift along the axial 
direction. 

There are a number of differences between turbu¬ 
lence in the wide-gap (e.g., radius ratio 0.5 to 0.8) and 
narrow-gap (e.g., radius ratio 0.99) Taylor-Couette sys¬ 
tems. Firstly, in the wide-gap case, the intensity distri¬ 
butions of turbulent fluctuations are often uneven in the 
radial direction 122 EZj, where more energetic turbulent 
fluctuations occur towards the inner cylinder wall. The 
regions near the inner cylinder thus exhibits stronger tur¬ 
bulence than in the region near the outer cylinder. This 
radial dependence of turbulent fluctuations is lost in the 
quasi two-dimensional, narrow-gap Taylor-Couette sys¬ 
tem, where turbulence is observed through the bulk in the 
radial direction. Secondly, in wide-gap systems, phenom¬ 
ena such as turbulent streaks [22] . small-scale Gortler 
vortices and herringbone-like streaks [23] can occur near 
both inner and outer cylindrical walls. Examining the 
typical size of the small Gortler vortices [23| reveals that, 
in the narrow-gap case these vortices have expansion 
larger than the radial width, excluding the possibility 
of generating turbulent streaks from small-scale Gortler 
vortices. Indeed, our simulations do not reveal any kind 
of such small scale vortices. This might explain the loss 
of radial dependence of turbulence, as can be seen, e.g., 
from the cross-flow energy behavior in Fig. [6] The mech¬ 
anism that generates turbulence in quasi two-dimensional 
Taylor-Couette system is thus quite different from that 
in the three-dimensional system. In fact, in narrow-gap 
systems turbulence emerges at the boundary layer of the 
neighboring vortex cells almost immediately at any radial 
position. 

Turbulence in quasi two-dimensional Taylor-Couette 
systems has different features than those in the planar 
systems as well. In particular, firstly, in planar Couette 
flows there is a significant difference between the veloc¬ 
ities at the two walls, and asymmetry in the intensity 


distribution of turbulent fluctuations is caused by the 
curvature effect. In the narrow-gap case, this difference 
is minimal. Secondly, while there are ring burst patterns 
of turbulent bands in planar Couette flows [HI? ], the 
background flow is laminar, versus wavy-like patterns in 
quasi two-dimensional systems, where the former defines 
a threshold between basic state and turbulence but the 
latter is a threshold between an already complex state 
and turbulence. Except this difference there is in fact a 
remarkable similarity between turbulence in both types 
of systems. For example, detailed investigations |3T of 
the laminar-turbulent boundary layer in planar Couette 
flows revealed isolated band states of turbulence in con¬ 
fined domains close to the global stability threshold. Al¬ 
though these bands appear at various different angles, 
they are all parallel, which are remarkably analogous to 
our ring-burst patterns on an unrolled cylindrical surface 
(e.g., comparing Fig. [l] with Figs. 1 and 11 in Ref. [? 
]). In addition, the routes to turbulence are similar: in 
both cases the turbulent bands grow out of a small local¬ 
ized turbulent spot that subsequently expands in some 
direction. 

It may be challenging to detect ring-burst patterns ex¬ 
perimentally as they coexist with other complex states 
such as VSWBs with similar turbulent characteristics. 
Nonetheless, given that the Taylor-Couette system is a 
paradigm enabling well controlled experiments on com¬ 
plex vortex dynamics and turbulence, we hope our find¬ 
ing will stimulate further research of turbulence in quasi 
two-dimensional regime of the system, a regime that has 
received little attention in spite of a large body of litera¬ 
ture on Taylor-Couette flows in general. 
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